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Thermal photon to dilepton ratio in high energy nuclear collisions 



o 
o 

(N 

m 



o 

=5 

(N 
> 

in 

7— I 

in 
o 

o 



X 



Jajati K. Nayak, Jane Alam, Sourav Sarkar and Bikash Sinha 
Variable Energy Cyclotron Centre, 1/AF, Bidhan Nagar , Kolkata - 700064 
(Dated: February 3, 2009) 

The ratio of transverse momentum distribution of thermal photons to dilepton has been evaluated. 
It is observed that this ratio reaches a plateau beyond a certain value of transverse momentum. We 
argue that this ratio can be used to estimate the initial temperature of the system by selecting 
the transverse momentum and invariance mass windows judiciously. It is demonstrated that if the 
radial flow is large then the plateau disappear and hence a deviation from the plateau can be used 
as an indicator of large radial flow. The sensitivity of the results on various input parameters has 
been studied. 

PACS numbers: 25.75.-q,25.75.Dw,24.85.+p 



I. INTRODUCTION 

Collisions between nuclei at ultra-relativistic en- 
ergies produce charged particles - either in the 
hadronic or in the partonic state, depending on the 
collision energy. Interaction of these charged parti- 
cles produce real and virtual photons (lepton pairs). 
Because of their nature of interaction, the mean free 
path of electromagnetic (EM) radiation is large com- 
pared to the size of the system formed after the col- 
lision. Therefore, EM radiation can be used as an 
efficient tool to understand the initial conditions of 
the system 0, HUSH 

and hence can be used to 
probe the quark gluon plasma (QGP) formation in 
heavy ion collisions (HIC). Practically, however, this 
is a difficult task, because on the one hand the ther- 
mal radiation from QGP has to be disentangled from 
those produced in initial hard collisions and from 
the decays of hadrons and on the other hand the 
evaluation of thermal photon and dilepton spectra 
need various inputs such as initial temperature (Tj), 
thermalization time (t^), equation of state (EOS), 
transition temperature (T c ), freeze-out temperature 
(Tf) etc, which are not known unambiguously. The 
sensitivity of the photon spectra on these inputs are 
demonstrated in Jy, 0] • Therefore, theoretical results 
on transverse momentum (pr) spectra of photons 
and dileptons always suffer from these uncertainties. 
Of course, certain constraints can be imposed on 
these inputs from experimental results - e.g. trans- 
verse mass spectra of hadrons and hadronic multi- 
plicities are useful quantities for constraining freeze- 
out conditions and initial entropy production. 

In the present work, therefore, we evaluate the 
ratio of the transverse momentum spectra of thermal 
photons and lepton pairs: 



Re 



{d 2 N 1 /d 2 p T dy) y=0 /(d 2 N 1 ,/d 2 p T dy) y=0 



in which most of the uncertainties mentioned above 
are expected to get canceled so that it provides ac- 
curate information [1, Q about the state of the mat- 
ter formed initially. We calculate the ratio, R em for 
SPS, RHIC and LHC energies. 

The paper is organized as follows. In section II 
we discuss the invariant yield of thermal photons 
and lepton pairs . In section III, the space time evo- 
lution is outlined. Results are presented in section 
IV. Finally in section V, we summarize the work. 



II. PRODUCTION OF THERMAL 
PHOTONS AND e+e PAIRS 

The rate of thermal dilepton production per unit 
space-time volume per unit four momentum volume 
is given by[l|, @, 0, % 
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a 



d 4 pd 4 x 12tt 4 M 2 



L(M 2 )ImU^f BE (2) 



a is EM coupling, Imll^ is the imaginary part of the 
retarded photon self energy and Jbe is the Bose- 
Einstcin factor which is a function of u^p^ for a ther- 
mal system having four velocity at each space- 
time point of the system, p 2 (~ p^p^) = M 2 is the 
invariant mass square of the lepton pair and 
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arises from the final state leptonic current involving 
Dirac spinors. m in Eq. [3]is the lepton mass. 

The real photon production rate can be obtained 
from the dilepton emission rate by replacing the 
product of the EM vertex 7* — > l + l~, the term in- 
volving final state leptonic current and the square of 
the (virtual) photon propagator by the polarization 



(1) sum (J2 



polarization 



-g^ u ) for the real photon. 
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Finally the phase space factor for the lepton-pairs 
should be replaced by that of the photon to obtain 
the photon emission rate, 
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d 4 xd 3 p (2tt) : 



rlmll^/BE 



(4) 



(see [JUSH for details )' 

The results given above is correct up to order 
e 2 (~ a) in EM interaction but exact, in principle, 
to all order in strong interaction. Now it is clear 
form Eqs.[2]and d]that for the evaluation of photon 
and dilepton production rates one needs to evaluate 
the imaginary part of the photon self energy. The 
Cutkosky rules or thermal cutting rules give a sys- 
tematic procedure to express the imaginary part of 
the photon self energy in terms of the physical am- 
plitude. 



A. Thermal photons 

Ideally, one wants to detect photons from QGP. 
However, the experimental measurements contain 
photons from various processes e.g, from the hard 
collisions of initial state partons of the colliding 
nuclei, thermal photons from quark matter and 
hadronic matter and photons from the hadronic de- 
cays after freeze-out. The contributions from the 
initial hard collisions of partons is under control via 
perturbative QCD (pQCD). The data from pp colli- 
sions will be very useful to validate pQCD calcula- 
tions. Photons from the hadronic decays (ir° — > 77, 
rj — > 77 etc.) can be reconstructed, in principle, by 
invariant mass analysis. But the most challenging 
task is to separate the thermal photons originating 
from the hadronic phase, which needs careful theo- 
retical estimation. 

The invariant yield of thermal photons can be 
written as 



d2pTdy ' 



d 

cPprdy 



d A x 



(5) 



where i = Q, M, H represents QGP, mixed (coexist- 
ing phase of QGP and hadrons) and hadronic phases 
respectively. (d 2 R/d 2 pTdy)i is the static rate of pho- 
ton production from the phase i, which is convoluted 
over the expansion dynamics through the integration 
over d 4 x. 



1. Thermal photons from Quark Gluon Plasma 

The contribution from QGP to the spectrum of 
thermal photons due to annihilation {qq^gj) and 



Compto n (q(q)fl q{q)l) processes has been calcu- 
lated in 10, 111 using hard thermal loop (HTL) ap- 
proximation [l2j. Later, it was shown that photons 
from the processes [13j : gq^gq^f, qq-^qq^, qqq^qj 
and gqq—>gj contribute in the same order 0(aa s ) as 
Compton and annihilation processes. The complete 
calculation of emission rate from QGP to order a s 
has been performed by resuming ladder diagrams in 
the effective theory [14j . In the present work this 
rate has been used. The temperature dependence of 
the strong coupling, a s has been taken from [Ta ] . 



2. Thermal photons from hadrons 

For the photon spectra from hadronic phase we 
consider an exhaustive set of hadronic reactions 
and the radiative decay of higher resonance states 
fl6l . Ir?! . Ojj]. The relevant reactions and decays 
for photon production are: (i) tt n — > p7, (ii) 
7r p — > 7T7 (with all possible mesons in the in- 
termediate state [lH) , (iii)7T7r — * 777 and (iv) 
ttt] — ► 777, p —>■ 77777 and uj — > 777. The corre- 
sponding vertices's are obtained from various phe- 
nomenological Lagrangians described in detail in 
Ref. [H, [ij, GJ|. The reactions involving strange 
mesons: ir K* — > K"f, ir K — > K* 7, pK — > A" 7 
and K K* — + 7r 7 [l!| have also been incorporated in 
the present work. Contributions from other decays, 
such as K* (892) Kj, <f> -> 777, &i(1235) -> 7r 7 , 
a2(1320) — > 777 and Ai(1270) — * 777 have been 
found to be small [H for p T > 1 GeV. All the 
isospin combinations for the above reactions and de- 
cays have properly been taken into account. The 
effects of hadronic form factors [l9[ have also been 
incorporated in the present calculation. 



B. Thermal dileptons 

Like photons, dileptons can also be used as an 
efficient probe for QGP diagnostics, provided one 
can subtract out contributions from Drell-Yan pro- 
cess, decays of vector mesons within the life time 
of the fire ball and hadronic decays occurring after 
the freeze-out. Like hard photons, lepton pairs from 
Drell-Yan processes can be estimated by pQCD. The 
Pt spectra of thermal lepton pair suffer from the 
problem of indistinguishability between QGP and 
hadronic sources unlike the usual invariant mass (M) 
spectra which shows characteristic resonance peaks 
in the low M region. The invariant transverse mo- 
mentum distribution of thermal dileptons (e + e~ or 
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virtual photons, 7*) is given by: 



d 2 N, 



7* 



d 2 p T dy 



E 

i=Q,M,H 



d 2 R^ 



d 2 p T dydM* 



dM 2 d 4 



(6) 

The limits for integration over M can be fixed ju- 
diciously to detect contributions from either quark 
matter or hadronic matter (see FigQ}. Experimen- 
tal measurements [2l], [22[ are available for different 
M window. 



the present work the space time evolution of the sys- 
tem has been studied using ideal relativistic hydro- 
dynamics in (2+1) dimension [27] with longitudinal 
boost invariance [28| and cylindrical symmetry. The 
initial temperature(Ti) and thermalization time (7$) 
are constrained by the following equation [2^] for 
an isentropic expansion: 
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1 dN 



45£(3) ia ef f ttR 2 a dy 
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1. Dileptons from QGP 

In the plasma phase the lowest order process pro- 
ducing lepton pair is qq — > 7* — > l + l~ . QCD cor- 
rections to this rate have been obtained for a QCD 
plasma at finite temperature in Refs. [23|, 0| up to 
order 0(a 2 a s ). In the present work contribution up 
to 0(a 2 a s ) has been considered. 



2. Dileptons from Hadrons 

The following parametrization [f| [2{| has been 
used to evaluate the dilepton emission rates from 
light vector mesons (p, lj and <p): 



d 2 IL 



dM 2 d 2 prdy 



a f \ 
1 



(M 2 - m 2 v ) 2 + (MT V ) 2 
1 



x(l 



8tt 1 + exp{{w - M)/8) 
a s 



-)]• 

7T 



(7) 



These parameterizations are consistent with the ex- 
perimental data from e + e~ — > V(p, uj or <fi) pro- 
cesses 0, H^l . Here, Jbe, is the Bose- Einstein 
distribution, fy is the coupling between the EM 
current and vector meson fields, my and IV are the 
masses and widths of the vector mesons and luq is 
the continuum threshold above which the asymp- 
totic freedom is restored. We have taken a s = 0.3, 
S = Q.2GeV, ujq — 1.3 GeV for p and u>. For <f> we 
have taken u Q = 1.5 GeV and 8 = 1.5 GeV. The EM 
current in terms of /?, to and <p field can be expressed 
as Jf, = JP + J£/3 - J^/3. Therefore, the contri- 
butions from uj and 4> will be down by a factor of 
9. 



where, dN/ dy= hadron multiplicity, Ra is the radius 
of the system, £(3) is the Riemann zeta function and 
a = Tr 2 g/90 (g = 2x8 + 7x2x2x3x N F /8) is the 
degeneracy of the massless quarks and gluons in the 
QGP, AV=mimber of flavours. The values of initial 
temperatures and thermalization times for various 
beam energies are shown in table I. The initial en- 
ergy density,e(xi, r) and radial velocity, u r (rj, r) pro- 
files are taken as: 



and 



v r (Ti,r) = vq I 1 - 



1 



(9) 



(10) 



where the surface thickness, 8 = 0.5 fm. We have 
taken vq — 0, which can reproduce the measured 
hadronic spectra at SPS and RHIC energies (5, l30l|. 
So far there is no consensus on the value of T c , it 
varies from 151 MeV [U to 192 MeV [32|. In the 
present work we assume T c = 192 MeV. In a first or- 
der phase transition scenario - we use the bag model 
EOS for the QGP phase and for the hadronic phase 
all the resonances with mass < 2.5 GeV have been 
considered [Uj]. 

To show the sensitivity of the results on the EOS 
we also use the lattice QCD EOS for T > T c [H. For 
the hadronic matter (below T c ) all the resonances 
with mass < 2.5 GeV have been considered [33(. For 
the transition region the following parametrization 
has been used [35j. 



a = f{T)s q + (1 - f(T))a h 



(11) 



where s q (s/j) is the entropy density of the quark 
(hadronic) phase at T c and 



III. SPACE-TIME EVOLUTION 



f(T) = \{l + tanh£^)) (12) 



The matter formed after ultra-relativistic heavy 
ion collisions undergo space-time evolution, which 
can be described by relativistic hydrodynamics. In 



the value of the parameter Y can be varied to make 
the transition strong or weak first order. Results for 
various values of V are given below. 
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TABLE I: The values of various parameters - thermaliza- 
tion time (ti), initial temperature (Ti), freeze-out tem- 
perature (Tf) and hadronic multiplicity dN/dy - used in 
the present calculations. 



Accelerator 


dN 

dy 


Ti{fm) 


Ti(GeV) Tf 


(MeV) 


SPS 


700 


1 


0.2 


120 


RHIC 


1100 


0.2 


0.4 


120 


LHC 


2100 


0.08 


0.7 


120 
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FIG. 1: The invariant mass distributions of thermal 
dileptons from QGP and hadronic matter at T — 200 
MeV. Solid (dashed) line indicates the emission rates 
from QGP (hadronic matter). The dot-dashed line 
stands for emission rate from hadronic matter at the 
transition temperature (see text). 



IV. RESULTS 
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FIG. 2: The thermal photon to dilepton ratio, Rem as a 
function of transverse momentum, pr for various invari- 
ant mass window. 
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The values of the initial and freeze-out parame- 
ters shown in table I along with the EOS mentioned 
above have been used as inputs to hydrodynamic 
calculations. The experimental data from SPS on 
hadrons 13611 . photons .37j and M distribution of 
dileptons |38| have been reproduced in [3(|, [39l.l4"oT| 
and (4l| respectively by using these inputs. The val- 
ues of the initial parameters for SPS agree with the 
results obtained from the analysis of photon spec- 
tra in Refs. 0, 53, 0, S| . Recently the data from 
PHENIX collaboration at RHIC [S| 53] has also 
been explained in Q (see also [48j |) with the parame- 
ters mentioned in table I. The lepton pairs measured 
by NA60 collaboration [Sil, [5(| in In-In collisions has 
been explained by spectral broadening of p [5l|, [52j ■ 

The emission rate from hadronic and quark mat- 
ter at a temperature of 200 MeV has been dis- 
played in Fig. [T] The contribution from QGP dom- 
inates over its hadronic counterpart (without any 
medium effects) for M < 600 MeV and M > 1.1 
GeV, therefore, these windows are better suited for 
the detection of QGP. However, it should be men- 



FIG. 3: Same as Fig. Q] for RHIC energy 
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FIG. 4: Same as Fig. Q] for LHC energy 
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FIG. 5: Same as Fig. [2] for quark matter phase only. 
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FIG. 6: Same as Fig. [2] for hadronic phase only. 



tioned here that the modification of the spectral 
functions of vector mesons (especially p and lu) - 
pole shift HU or broadening [26| may give rise to 
dileptons at the lower M region making it difficult to 
detect contributions from QGP. The change in the 
hadronic spectral function will enhance the dilep- 
tons from the hadronic contribution in the lower 
mass (M < 600MeV) window, however the over- 
all structure in the ratio, R em will not change ap- 
preciably. At the transition temperature (~ 200 
MeV) if one assumes the vector mesons masses go 
zero a la Brown-Rho scaling [53] then all the peaks 
in the dilepton spectra disappeared and the rates 
obtained from EM current-current correlator (dot- 
dashed line) are close to the rate from QGP, indi- 
cating that the qq interaction in the vector channel 
has become very weak, signaling the onset of de- 
confinement. This also indicates the quark- hadron 
duality 0, near the transition point. 



FIG. 7: The variation R em with pr for invariant mass 
window, M — 0.7 — 0.8 GeV. An unrealistically large 
value to radial flow has been given initially to demon- 
strate that large flow can destroy the plateau structure 
of R e m- Other inputs are similar to those of Figs[3]and 
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FIG. 8: Rem as a function of pr for different values of 
T c for invariant mass windows, M — 0.7 — 0.8 GeV and 
M = 1.2 - 1.3 GeV. 



It is well known that the px spectra of photons and 
lepton pairs are sensitive to the values of initial tem- 
perature Ti, vq, Tf and EOS. The T c dependence of 
the pt distribution is found to be negligibly small [|| . 
As we have mentioned before though these parame- 
ters can be constrained from the measured multiplic- 
ity and freeze-out spectra there remains still some 
room to vary these quantities in order to be able to 
describe the experimental data. 

The pt dependence of the ratio, R em for SPS, 
RHIC and LHC energies are shown in Figs. [2 [3] and 
[3] respectively. It is observed that at a given px, the 
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FIG. 9: Rem as a function of pt for different EOS for 
invariant mass windows, M = 0.7 — 0.8 GeV and M = 
1.2 - 1.3 GeV. 
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It is clear from the results displayed in Figs. [H [3] 
and [4] that the quantity, R em , reaches a plateau 
beyond pr = 1-5 GeV for all the three cases i.e. for 
SPS, RHIC and LHC. It may be noted here that the 
degree of flatness increases from SPS to RHIC and 
LHC. As mentioned before for all the three cases, 
except Ti all other quantities e.g. T c , vq and EOS 
are same, so the difference in the value of R em in the 
plateau region originates due to different values of 
initial temperature, indicating this can be a measure 
of Ti. 

The following analysis will be useful to under- 
stand the origin of the plateau at high px re- 
gion. The strong three momentum dependence in 
the dilepton and photon emission rates (Eqs. and 
d] respectively) originates from the thermal factor, 
fBE(E, T). For a static system the energy, E can be 
written as E = Mrcoshy, where Mt = \Jv\ + A/ 2 
y = tanh- 1 p z /E. At high p T (» M), M T » 
Pt, the exponential momentum dependence become 
same for real photon (M 2 = 0) and dilepton (M 2 ^ 
0) spectra and hence plateau is expected in the static 
ratio for large pr for all the M values. 

We recall that for an expanding system out of the 
two kinematic variables describing the dilepton spec- 
tra, pt is affected by expansion but M remains un- 
changed. The range of M under present study is 
0.3 < M(GeV) < 1.3. The energy, E appearing in 
both the photon and dilepton emission rates should 
be replaced by u^p^ for a system expanding with 
space-time dependent four velocity u M . Under the 
assumption of cylindrical symmetry and longitudi- 
nal boost invariance can be written as 



u 1 * = "f r (t/T,v r cos(j>,v r sin(j),z/T) 



(13) 



FIG. 10: Initial temperature is plotted as function 
Re m {pT = 2.5GeV) for the M window, M = 1.2 - 1.3 
GeV. 



where r = \/t 2 — z 2 t = rcoshrj, z 
rsinhrj, v r (r,r) is the radial velocity, 7 r (r, r) 



(1 — v r (T,r)) 1 I 2 . The four momentum, 



(MTCoshy,pT,0, MTsinhy) where pl 
Therefore, for dilepton 



rriTsinhy. 



ratio decreases with M, reaches a minimum around 
p-peak and increases beyond the p-peak. This trend 
is valid for all the cases, i.e. SPS, RHIC and LHC as 
expected because at a given pt the R em is actually 
the inverse of the invariant mass distribution of lep- 
ton pairs (the denominator i.e. the photon spectra is 
same for all the mass windows). It is observed from 
Figs. and Q] that the ratio, R em decreases with 
Ti for given px for M below the p-peak and the op- 
posite behaviour is observed above the p-peak. The 
slope of the ratio at low pt also indicates substan- 
tial change with increasing M, the slope is minimum 
at the p-peak. Therefore, the minimum of the slope 
may be used to locate the effective mass of the vector 
meson in medium. 



u^Pv 



j r (MTCOsh(y — 77) — v r pTCos4>) (14) 



for photon the factor u^p^ can be obtained by re- 
placing Mt in Eq. [14] by px- The pt dependence 
of the photon and dilepton spectra originating from 
an expanding system is predominantly determined 
by the thermal factor /be- Therefore, we discuss 
following three scenarios. (i)At high pt{» M), 
Mt ~ Pt, the exponential momentum dependence 
become same for real photon and dilepton spectra, 
hence for large pt a plateau is obtained in the ratio, 
Rem (Figs [2HU ■ In other words, the effect of radial 
flow on the photon and dilepton is similar at high pt 
region. (ii)If the large M pairs originate from early 
time (when the flow is small) the ratio, R em which 
includes space-time dynamics will be close to the 
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static case and hence will show plateau, (iii) How- 
ever, at late time when the radial flow is large and 
M is comparable to or larger than pt the effect of 
flow on dilepton will be larger (receives larger radial 
kick due to non-zero M) than the photon and hence 
the plateau may disappear. Therefore, the disap- 
pearance of plateau structure in R em in moderate 
or high M region will indicate the presence large ra- 
dial flow. This can be understood from the results 
shown in Figs [5] and [6] 

In Fig. [5] the ratio has been displayed only for 
quark matter. Here the flow is expected to be lower 
within the present framework of present study. A 
plateau is observed for all the M windows. It is 
observed that for high M (~ 1.2 GeV) and low M 
(~ 0.3 GeV) the ratio for QM is close to the total 
for LHC energy (not shown separately). 

In Fig[6]the ratios has been displayed for hadronic 
matter only. Here the flow is expected to be very 
large. Within the ambit of the present modeling 
the contribution from the hadronic matter is over- 
whelmingly large in the M region, 0.7 < M < 0.8 
GeV. Therefore, this region will have large effects 
from the radial flow and hence it may destroy the 
plateau. This is clearly seen in Fig. [6] for the curve 
corresponding to 0.7 < M < 0.8 GeV. 

To demonstrate the effect of flow on the plateau 
we use an initial velocity profile (which gives rise 
to stronger radial flow than Eq. [T0|) of the form 
v r (Ti,r) = v' jf— with an unrealistically large value 
of v' Q ~ 0.5. These inputs are used only for results 
shown in FigJ7J which clearly indicates the disap- 
pearance of plateau. Variation of R em with px cor- 
responding to hadronic phase is steeper than the to- 
tal because of larger radial flow in the late stage of 
the evolution. 

Now we demonstrate the effect of other parame- 
ters on Rem- The value of T c has large uncertainties. 
Therefore, we show the sensitivity of the results on 
T c in Fig. [5] for two invariant mass windows. The 
results are insensitive to T c . 

The effect of the EOS on R em is demonstrated in 
Fig. O by varying T in Eq. [T3J It is observed that 
the effect of EOS on R em for both the mass windows 
are small. Similar to the effect of T c , the larger mass 
window (1.2 < M(GeV)< 1.3) is less affected by the 
change in EOS. This is because the effect radial flow 
(and other hydrodynamic effects) are less at early 
times from where higher mass lepton pairs originate. 
Replacement of lattice QCD EOS for QGP phase by 
bag model shows negligible effects on R em . 

In Fig.rrU]the dependence of R em (pT — 2.5GeV) is 
depicted as a function of T t for 1.2 < M(GeV)< 1.3. 
This mass window is selected because the contribu- 
tions from the hot quark matter phase dominates 
this region and the effects of T c , EOS etc are least 
here, pr = 2.5 GeV is taken because R em achieved 




FIG. 11: The variation R em for hard photons to dilep- 
tons ratio as a function of pr for ^/sjvjv = 200 GeV and 
invariant mass window, M = 0.2 — 0.3 GeV. 



a complete plateau at this value of transverse mo- 
mentum. The change in R em from SPS to RHIC 
is about 40% and from RHIC to LHC this is about 
20%. A simultaneous measurements of photons and 
dileptons with required accuracy, will be useful to 
disentangle the effects of flow and true average tem- 
perature in a space-time evolving system formed in 
heavy ion collisions at ultra-relativistic energies. 

We have evaluated RP m CD , the ratio 
(d 2 N 1 /d 2 pTdy) y= o/(d 2 N~ ( */d 2 pTdy)y = o for hard 
processes using pQCD (Fig. Hip. The hard photon 
contributions has been constrained to reproduce the 
PHENIX data [1| for pp collisions at ^/sjf^ = 200 
GeV. We consider qq — > 7* — > l + l~ 7 qq^gj* and 
qg(q)—>qq"f* for the lepton pair production. The M 
integration of lepton pair spectra (Eq. [S]) is done 
over the range 0.2 < M(GeV) < 0.3. We observe 
that R v em CD increases for pr up to ~ 3 GeV, above 
which it reaches a plateau. Therefore, for px ~ 1 — 3 
GeV, Rem for the thermal and pQCD processes 
show different kind of behaviour. The plateau arises 
from the fact that at large px both photon and 
dilepton show power law behaviour [5 71 1 5 81 ] . In the 
low pt domain lepton pairs (photon) from pQCD 
processes indicate a Gaussian type [57| (power law) 
variation resulting in the increase of R em with px- 

V. SUMMARY AND CONCLUSIONS 

We have studied the variation of R em , the ratio 
of the transverse momentum spectra of photons and 
dileptons and argued that measurement of this quan- 
tity will be very useful to determine the value of the 
initial temperature of the system formed after heavy 
ion collisions. We have observed that R em reaches a 
plateau beyond pt = 1.5 GeV. The value of R em in 
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the plateau region depends on T^. However, the ef- 
fects of flow, EOS and the dependence on the values 
of T c , vq and other model dependences get canceled 
away in the ratio, R em - For M above and below the 
p peak and px > 2 GeV the contributions from quark 
matter dominates, therefore these regions could be 
chosen to estimate the initial temperature of the sys- 
tem formed after the collisions. 

It is well known that T e fj, the inverse slope (see 
e.g. extracted from the pr spectra of EM ra- 
diation contains the effect of temperature as well as 
flow. We have seen that when the flow is less (in the 
initial stage of the evolution) the ratio, R em shows 
a plateau for large pr{» M), the height of the 
plateau in this region will give a good measure of 
the average temperature. However, a large flow can 
destroy the plateau and hence the deviation from 
the flatness of the R em versus pr curve may be used 
as a measure of flow. So a careful selection of M 
and pt regions will be very helpful to disentangle 
the effect of true average temperature and the flow 
(see also [5^]). In (5t| it was shown that the effects 
of flow on high M(> 1.5 GeV) could be quite large 
and in such cases the plateau in R em may disappear. 
However, in the present work we confine in the range 
0.3 < M)GeV) < 1.2. 

EM radiations originating from the interactions 
between thermal and non-thermal (high energy) par- 
tons [gO, [HI, has been neglected in the present 
work. It is expected that the EM radiation from 
these processes and also from the pre-equilibrium 
stage will not affect R em . 

We have studied the effects of chemical off- 
equilibrium of mesons on the photon and dilepton 
production rates. This is implemented by appro- 
priately introducing non-zero pionic chemical po- 
tential, p n (p p = 2/1^,/^ = 3/i7r) in the thermal 
factors [63| appearing both in photon and dilepton 
emission rates. We observed that the plateau struc- 
tures in R em do not change for RHIC and LHC, but 
for SPS it has little effect. 



The change in hadronic spectral function at non- 
zero temperature and density is a field of high con- 
temporary research interest as this is connected with 
the restoration of chiral symmetry in QCD. From 
the QGP diagnostics point of view the background 
contributions (photons and dilcptons from thermal- 
ized hadrons) are affected due to medium effects on 
hadrons. Therefore, some comments on this issue 
are in order here. 

We have checked that the pr spectra of both pho- 
tons and dileptons are sensitive to the pole shift 
of hadronic spectral function, as the reduction of 
hadronic masses [53[ in a thermal bath increases 
their abundances and hence the rate of emission gets 
enhanced 0, [lfl [I?], H^, H(| • The invariant mass dis- 
tribution of lepton pairs are sensitive to both the 
pole shift and broadening [H, [H EH, [H, HI . But 
the pt spectra of the EM radiation is insensitive 
to the broadening of the spectral function provided 
the integration over the M is performed over the 
entire region. This is because broadening does not 
chan ge t he density of vector mesons significantly (see 
also [40j). However, the number density of vector 
mesons depends on the nature (shape) of the spec- 
tral function within the integration limit. There- 
fore, the pt spectra may change due to broadening 
when the integration over M is done in a limited 
M domain. We have checked that doubling the p 
width (~2x 150 MeV) changes R em by 10%. It 
is important to note that the change in mass and 
widths can not be arbitrary it should obey certain 
constraints as discussed in [66]. Therefore, simul- 
taneous measurements of pt spectra and invariant 
mass distribution of real and virtual photons could 
be very useful to understand the nature of medium 
effects on hadrons [40l |. 
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